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Fluctuating entropy production is studied for a set of linearly coupled complex fields. The general 
result is applied to non-equilibrium fluctuating hydrodynamic equations for coarse-grained fields 
(density, temperature and velocity), in the framework of model granular fluids. We find that the 
average entropy production, obtained from the microscopic stochastic description, can be expressed 
in terms of macroscopic quantities, in analogy with linear non-equilibrium thermodynamics. We 
consider the specific cases of driven granular fluids with two different kinds of thermostat and 
the homogeneous cooling regime. In all cases, the average entropy production turns out to be 
the product of a thermodynamic force and a current: the former depends on the specific energy 
injection mechanism, the latter takes always the form of a static correlation between fluctuations 
of density and temperature time-derivative. Both vanish in the elastic limit. The behavior of the 
entropy production is studied at different length scales and the qualitative differences arising for the 
different granular models are discussed. 
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I. INTRODUCTION 



Among the several different efforts to describe out-of-equilibrium systems, the linear non-equilibrium thermody- 
namic approach [l| has obtained many important results, at least for a class of systems close to equilibrium. In this 
framework, a fundamental role is played by the macroscopic entropy production, which is related to the irreversible 
phenomena occurring in the system due to the presence of heat and matter currents, chemical reactions, viscous flows, 
etc... For systems in contact with a single thermal bath, the study of fluctuations relaxing towards the equilibrium 
state can be described within such a theory, and extensions of it 0. Moreover, the great interest in the study of 
entropy production in non-equilibrium stationary states has been motivated by the fact that, in some particular cases, 
it can be shown that the stationary state is characterized by the minimum entropy production principle 0, Q . 

For generally far from equilibrium systems, a comprehensive theory is still lacking, but some important results 
have been obtained in the last years. In particular, in the framework of systems described by stochastic models, a 
stochastic thermodynamics has been proposed, extendi ng th e concepts of standard thermodynamics to fluctuating 
quantities, and some general relations have been proved [4l4l3l|. As for classical equilibrium thermodynamics, a central 
role is played by the microscopic, fluctuating version of entropy production, defined as functional of finite time-length 
trajectories in the phase space. Such a quantity satisfies relations which have been considered by some authors as 
generalizations of the second law of thermodynamics [l4| • 

The specific aim of this paper is the study of fluctuations of entropy production in systems described by a set of 
coupled Langevin equations. Many examples of physical systems are well characterized by such a kind of stochastic 
process [15l - ll8l |. In particular, we are interested in non-equilibrium stationary systems where the coupling among the 
different fields defining the model, as well as the noise entering the Langevin description, violate the detailed balance 
condition. In these out-of-equilibrium systems a finite rate of entropy production can be measured. However, from a 
microscopic point of view, fluctuating entropy production is an elusive quantity: in a non-equilibrium stationary state 
(NESS) the total entropy S(t) (assuming we have a definition for it) has vanishing time-derivative. In analogy with 
the perspective of linear non-equilibrium thermodynamics, one can try to separate this zero time-derivative into two 
opposite contributions: one is a surface term, namely the entropy flow through the boundaries (e.g. heat exchanged 
with the thermostat) and the other is a bulk, volume term, namely the entropy produced inside the system because 
of internal lack of detailed balance: 



dS__ dS_ 
~dt ~ ~dt 



dS_ 
~dt 



(1) 

prod 

At equilibrium both terms are vanishing, but in a NESS the production term is always positive. In linear non- 
equilibrium thermodynamics, symmetry arguments contribute to make this separation unique; however, in general, it 
is not evident that this unique separation can always be achieved. A definition of the entropy production term related 
to the dynamics rather than to the thermodynamics of the system has been proposed in recent years: 

dS 



dt 



= MM> o, (2) 

, t— >oo t 
prod 



where E(£) is a functional of the trajectory {4>(s)} of the system in the time interval s £ [0, t], with 4>(s) a set of relevant 
observables. As will be also shown in Sec. IIV( (£) is often associated with particular currents (or time-asymmetric 
correlations) vanishing at equilibrium. Fluctuations of S on finite-time trajectories in finite-size systems, may display 
also negative values, with an exponentially small probability, as dictated by the Fluctuation Relations @, 

In this paper we approach the problem on fairly general grounds. Within the framework of Markovian dynamics, we 
bring to the fore the structural elements which are common to general systems of Langevin equations. In particular, 
we relate the expressions obtained from the stochastic description of the system to macroscopic quantities, which can 
be put in a form very similar to that of linear non-equilibrium thermodynamics. 

In order to illustrate our gen eral results, we shall focus on the set of non-equilibrium fluctuating hydrodynamic 
equations for granular fluids [l9l[20j. We consider both driven granular gases, where two different kinds of thermostat 
are studied, and the case of homogeneous cooling (2lj . For such systems, in certain ranges of the physical parameters, 
a description in term of coarse-grained hydrodynamic fields can be given, and a set of coupled Langevin equations can 
be written, involving density, temperature and velocity fields. Applying our formalism to these systems, we are able 
to express the average entropy production in terms of macroscopic currents and thermodynamical forces. A similar 
result has been recently obtained in [221 ] . in the complementary framework of non-equilibrium stationary spin systems, 
described by transition rates violating detailed balance, and in [l2l |. 

The structure of the paper is the following. In Sec. [TT] we briefly recall some concepts of linear non-equilibrium 
thermodynamics. In Sec. IIIII wc obtain an expression for the entropy production of a set of Langevin equations, which 
is applied to fluctuating hydrodynamic equations for granular fluids in Sec. IIV1 Finally, some conclusions are drawn 
in Sec. [V] Two Appendices are devoted to some technical details. 
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II. ENTROPY PRODUCTION IN LINEAR NON-EQUILIBRIUM THERMODYNAMICS 

Here we recall some concepts of the linear non-equilibrium thermodynamics theory [l[ , in order to give a reference 
frame for the following discussions. 

Linear non-equilibrium thermodynamics is a continuum theory which aims at describing macroscopic systems 
characterized by irreversible processes. The starting point is the balance equation for the entropy, stating that the 
entropy in a volume changes because entropy flows from the boundaries into the volume, or because some irreversible 
phenomena are taking place inside the volume. The underlying hypothesis is a "local" equilibrium assumption, which 
allows one to write the thermodynamic Gibbs relations connecting the local (within a small mass element) entropy 
with other thermodynamic quantities. 

A. Thermodynamic forces and fluxes 

The first step of the theory is the local formulation of the conservation laws. In general, the conservation law for a 
given conserved quantity p (energy, mass, a component of momentum) in the system reads 

where J p is the flux associated with p. For an open system, in the presence of sources (or sinks) for the quantity p, 
we write the balance equation 

9 pit) 

where v p represents the production (or absorption) of quantity p in unit time. 

As mentioned above, the variation of entropy S of a macroscopic system is described by a balance equation, namely 
it can be split in two contributions: 

dS=dS\ ext +dS\ prod , (5) 

where dS\ ex t is the entropy change due to the coupling of the system with the surrounding medium, and dS\ pro d is 
the entropy produced inside the system due to irreversible processes. For an insulated system, dS\ ex t = and then 
dS = dS\ pro d > 0, from the second law of thermodynamics. For closed systems exchanging heat Q with a thermostat 
at temperature T, one has dS\ ext = ^ and dS > ^r. The change in time of the total entropy S is 



dS__ dS_ 
~dt ~ ~dt 



dS 
~dt 



J s • dA + / s dV, (6) 

prod J A JV 



where A denotes the contour surface of the system, J s is the entropy flow for unit time and unit area, s is the entropy 
production density per unit time and V the total volume. For open systems, J s is the sum of the heat flow divided 
by the local temperature, J q /T, plus all the flows of matters from the outside. Using the conservation laws and 
"local" equilibrium hypothesis, the entropy production can be related to the several different irreversible phenomena 
occurring inside the system. The structure of the entropy production s is then that of a bilinear form: 

* = £>-Fi, (7) 

i 

where Ji are fluxes (or currents) of the quantities describing the system, which are associated with irreversible 
phenomena, and Fi are thermodynamic forces (or affinities), related to gradients of state variables (spatial non- 
homogeneities) or to external forces. In equilibrium conditions, all the thermodynamic forces vanish and so the 
entropy production does. Moreover, one also requires that all fluxes vanish with the thermodynamic forces. 

Let us consider for instance the case of a metal bar coupled at the edges with two thermostats at different temper- 
atures [f2[. In this situation, a stationary temperature profile sets up along the bar, with a surface entropy flux due 
to flux of heat across the edges and a bulk entropy production due to the sustain of a temperature gradient in the 
bulk: 

1 VT 
,s=-J q -F q = -J q -^ r , (8) 
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where T is the local temperature and the thermodynamic force conjugated to the heat flux is F q = — (VT)/T 2 . 
From the entropy production in the bulk and the knowledge of the geometry of our problem the net flux of entropy 
across the edges can be then deduced. In Sec. IIV1 we shall consider homogeneously driven granular systems in two 
dimensions, where such a simple distinction between surface and bulk contributions cannot be carried out. 



B. Onsager coefficients 



In order to close the system of equations for the entropy production, conservation laws and entropy balance equations 
have to be supplemented by the phenomenological relations, which, as first approximation, express the fluxes in terms 
of a linear combination of thermodynamic forces 



through the Onsager coefficients Ly. Substituting such relations into the expression for entropy production, one 
obtains 



* = £ (10) 



namely a quadratic expression in the thermodynamic forces. In equilibrium conditions, the thermodynamic forces 
vanish and, due to Eq. ©, so the fluxes do. 

In the physical example introduced above, the phenomenological relation is the Fourier's law 



VT 

Jq = -L qq j^ = — AVT, (11) 



where A is the heat conductivity, and 

V T 2 J 



s = L qq (?Pl . (12) 



If the forces Fi perturb a system at equilibrium, i.e. a system where only conservative (and therefore time-reversible) 
forces are present, then Onsager has shown that the coefficients satisfy the so-called reciprocal relations {231 l24j) . 



III. FLUCTUATING ENTROPY PRODUCTION IN STOCHASTIC DYNAMICS 



Statistical systems are conveniently described by stochastic processes. We recall here how the notion of entropy 
production can be introduced in such a framework. In particular, here we focus on the degree of irreversibility of 
the dynamics, as measured by the functional S(^o) introduced by Lebowitz and Spohn Q, where fig = {<f>(s)} is a 
trajectory of the phase space in the time interval s 6 [0,i]. Let us consider for simplicity a discrete time process, 
where jumps occur at times U, with i € [0,n]. The functional S(f2 ) is defined as the ratio between the probability 
of a given trajectory P(f2 ) = p(<l>o)W (<fro\(/>i) ■ • ■ W(<fr n -i\<frn), where (pi = <p(ti), and the probability of time-reversed 
trajectory T(£l ): 



£(t) = ml) 



In 



p(</>o) 



In 



W{M<l>i)---W{<P n -i\cp n ) 



(13) 



where f2 



{4>{t — s)} and <f> = e(j), with e = ±1, according to the parity of the field. 
In a stationary state, where entropy produced/consumed inside the system is continuously balanced by a flux 
coming from the boundaries, as in Eq. ([6]) 







dS 
~dt 



dS_ 
~dt 



dS_ 
~dt 



(14) 



prod 



we write, following Q, 
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dS_ 
~dt 



= lim - (E(t)) > 0. (15) 

, t— foo t 
proa 



The functional E(t) depends on a stochastic variable (i.e. a trajectory of a stochastic differential equation) and is 
itself a stochastic variable with a given probability density: indeed, for finite t, it can also take negative values, of 
course with a probability decreasing with increasing t. At large times and for bounded (or negligible) values of the 

term In , the probability density of E satisfies the so-called Fluctuation Relation 0, H, 0, EH [H, [26| . 

In the following we will show that, in cases which are relevant for non-equilibrium spatially extended systems, the 
functional (|13[) can take a clear thermodynamic meaning, since it can be expressed in terms of macroscopic currents 
and forces, in strict analogy with the form of entropy production (J7)) discussed in Sec. Ill Al 



A. Coupled Langevin equations 



Here we specialize to a continuous space and time Markov process and, hereafter, </> denotes a complex d$- 
dimensional vector, with components fa (i = 1,...,^). The time-reversal transformation on vector <f> is defined 
as an operator which changes fa — > fa = eifa with G {+1,-1}: this implies that both real and imaginary parts of 
fa have the same sign-change upon time-reversal. As shorthand notations we will use <f> or e(f> to indicate the vector 
made of time-reversed components {eifa}. In a coarse-grained description of a fluid system, instances of even variables 
are the Fourier components of local density and temperature fluctuations (q = +1), whereas the velocity field is odd 
(f j = —1). We assume that the dynamics of each component of the vector </> is described by a Langevin equation: 



^-= Bi [fat)}+m\ (i6) 

with Bi a complex drift and & a Gaussian process with (£i(£)£J (s)) = 2Sij6(t — s)Da, where the x* denotes the 
complex conjugate of x. If we consider a trajectory between time and time t, we can define its time-reversal as 
fa(s) = tifait — s). For the complementary discussion starting from the corresponding Fokker-Planck equation, see 
Appendix [Bj The stochastic system in Eq. (jTHJ) is indeed a good model to describe the behavior of macroscopic 
variables in many contexts fl6j . Often it works as a first approximation, where one retains only the linear part of the 
dynamics and the effect of nonlinearitics is replaced with noise terms. 

In order to define a useful projection of the dynamics, let us for the moment neglect the noise terms in Eq. (|16[) . 
namely = 0. The time-reversal trajectory satisfies the following differential equation: 



dfajs) _ _ dfa (t - s) 
ds dt 

We notice two particular cases for Bi(fa): 



-e l B l [fat-s)] = -e i B i [efa[s)]. (17) 



B % {efa = -tiBiifa) (18) 
B t {efa = e^(0). (19) 



In the first case fa(s) satisfies exactly the forward equation (fT6|) . In the second case it satisfies the same equation 
with drift changed of sign. Following these two limit cases, we will in general decompose the i-th component of the 
drift as [HI 

Bi{fa) = B i>rev (<t>) + B i>irr (cf>), (20) 

where 

Bi, rev {4>) = \[Bi{fa) - eiBi{efa)] = -e.B^efa (21) 
Bi,irr{fa) = -x[Bi{fa) + eiBiiefa)] = €iB iiirr (efa, (22) 
and Bi t i rr represent the reversible and irreversible contributions to the drift, respectively. 
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B. Entropy production functional 

We assume that the system (| 16[) converges to a unique stationary state in the long time limit. Then, let us consider 
the entropy production in such a stationary state: following the Onsager-Machlup prescription, in the case of an 
additive Gaussian noise wc can write the probability of a trajectory J7q = {(f>(s), s = ■ • • t} as: 



P(nl) = p[#0)] exp[5(n* )], P(n*) = p[#0)] «*p[S(n$)], (23) 
where p(<f>) is the weight (probability or density) of state <j> in the steady state, and the action S is given by [27j 



5(Q* ) 



Jo 



AD. 



ii,R 



d(f>j{s) 
ds 



BMs)] 



ds 



b; [</>(*)] 



(24) 



We have used the notation Y^'i to indicate that the sum runs only on those indexes such that Da ^ [28]. Due 
to stationarity, we can always set to zero the initial time of the interval where the entropy production is calculated. 
Recalling that J* dsf(t — s) = f Q dsf(s), we easily get for the entropy produced by a path SIq [29t | : 



ml) = In - In EM = S(nl) - S(&) 

p(ng) i#(t)] 



(is cr(s), 



(25) 



with entropy production rate 



a(s) 



Bi,irrW)\ (j>*M-Bi re Ms)])+C.< 



(26) 



force 



where c.c. denotes the complex conjugate. In this formalism, we can point out how the entropy production is expressed 
as the product of forces (the drift) multiplied by currents (time derivative of fields). According to the phcnomcnological 
relations (|9]), the irreversible fluxes vanish when the driving forces are switched off. Our identification of forces and 
fluxes is consistent with such relations: indeed, when (B iyirr [(j)(s)}) = we also have (</>*(s) — B* rev [cf)(s)]) = 0, as 
can be seen by averaging equation (|16|) over the noise. This is in analogy with the macroscopic expressions in linear 
non-equilibrium thermodynamics. 



We consider now the linear case B[4>] 



C. Linear processes 

where A is the so-called dynamical matrix: 



B i [<j>]=Y,A i j<l>i= E A v&+ E 'Mi' ( 27 ) 

j jGrcv(i) jeirr(i) 

where the shorthand notation j £ rev(i) means that the index j runs on the set of indices such that tj = — Cj, while 
j e irr(i) stands for the set of indices such that Cj — £j. In this case one has 



- = E^-f E 



.4 



«7 VO 



i jeirr(i) 



<t>*- E A *i 



j'Grev(i) 



+C.C. 



(28) 



force 



t- 

^ 2D,, 



jeirr(i) 



Jgirr(i) 



(29) 



In Appendix |A1 we show that at equilibrium (a) = 0, as expected. Formula (|28[) is analogous to the macroscopic 
thermodynamic result ([7]): it expresses the entropy production as the product of a force by a current. However, in 
the general case, the formula remains rather abstract, and we need explicit examples to illustrate its meaning. These 
will be discussed in the following section. 
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IV. APPLICATIONS TO GRANULAR FLUIDS 

In this section we show three examples in the framework of granular fluids where, starting from the definition (|25[) . 
we can obtain an expression for the entropy production similar to Eq. ([7]). 

A. Fluctuating hydrodynamics for driven granular fluids 

We consider a fluid of N identical inelastic hard spheres in dimension d, of mass 1 and diameter r, in a square 
box of volume V with external homogeneous stochastic driving. We denote by p the packing fraction of the system, 
which for d = 2 reads p = Nir(r /2) 2 /V . The model is defined by giving the equation of motion for the velocity of 
i-th particle: 



Vi (t) = -7bVi (t) + VZTWbQ (t) + Fi(t), (30) 

where 7& is a viscous drag, T& is the temperature of the external thermostat, Q is a Gaussian noise with (£;) = 
and {(i{t)Q(t')) = Sij5(t — t') and Fi(t) is the resulting force of eventual instantaneous collisions with other particles. 
Every time two particles i and j collide, their velocities are instantaneously changed following the rule 

1 + a 

v t =Vi — ~ v 3 ) ■ n]h, (31) 

where a G [0, 1] is the restitution coefficient (= 1 for clastic collisions) and h is the unit vector in the direction joining 
the centers of the colliding particles. This model has been introduced in [3(| and, recently, it has been also successfully 
used to describe real granular experiments [3l|, [Hj]. A different model thermostat is obtained in the limit ■fb — s- 0, 
Xh — > oo with constant T = 2^bTb, as it has been studied in [l9[ (see Sec. IIV A3p . In both cases the fluid reaches a 
NESS with a well defined granular temperature T = (\v\ 2 )/d. In the case of finite jb > 0, a characteristic thermostat 
time Tb = I/76 is defined. When compared to the mean collision time r c = (with to the collision frequency), it 
determines two possible regimes: Tb 3> t c is the dissipative regime, where T < Tb (the equal sign holding only for 
a = 1); T(, <C t c is the equilibrium regime where T = Tb for any a, since collisions are negligible. 

Let us add a brief comment on such a model, slightly anticipating some considerations on the entropy production 
to be discussed below. From Eq. (f3"0)) . we see that the thermostat is homogeneously coupled to all the particles. 
Therefore, the energy injection represents by definition a bulk contribution. It is then through collisions, Eq. (|3ip . 
that the dissipation of energy takes place. Because collisions occur across the whole system, also the energy sink 
represent a bulk contribution. Therefore, we can imagine that inelastic collisions are perfectly equivalent to a zero 
temperature reservoir, which is coupled to all the particles. In such a description, what is generally called "entropy 
flux" refers to the entropy exchanged with reservoirs, rather than to entropy flowing into the system through the 
boundaries. 

In the case where a separation of temporal and spatial scales takes place, a hydrodynamic description of the 
system can be given. Then, we introduce the coarse-grained hydrodynamic fields, density, velocity and temperature, 
n(r, t), u(r, t) and T(r,t), respectively, where r denotes a point in the d-dimcnsional space, as follows: 

n(r,t)=53<y(r-r i (t)) I 

i 

u{T,t) = -Y i v i (t)S(T-T i {t)), (32) 
n z — 4 

i 

i 

Here the sums are over all the particles in the system. The homogeneous stationary state is characterized by constant 
density n and granular temperature T. The fluctuations of hydrodynamic fields are {Sn = n(r, t) — n, ST = T(r, t) — 
T, u x , u y } and are described by the set of linear hydrodynamic equations. In these equations the small scale fluctuations 
have been projected out, but their effect on large scale fluctuations can be taken into account by a proper addition 
of noise terms UM, resulting in a fluctuating hydrodynamic description. For these models, such equations have been 
studied in fl9l |2G|. 

After linearization near the homogeneous state, and changing to Fourier space, in two spatial dimensions (d = 2), 
we are left with a linear Langevin system for the components of a four-dimensional (d^ = 4) complex field vector <j) = 
(5n(k), 5T(k), u{k), v(k)), with parities under time-reversal e — (1, 1, — 1, —1). By u and v we mean the longitudinal 
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and transverse velocity field, respectively. In this case a useful shorthand notation consists in replacing the index i 
by a label equal to the name of the fields (omitting the 5s for simplicity), i.e. i = 1 — > n, i = 2 — >T, i = 3 — > u and 
i = 4 — > v. The Langevin equation for each component of the vector field is then 



(33) 



where the specific form of the dynamic matrix A and of the noise amplitudes depend on the kind of thermostat and 
will be explicitly given below. In this case, from the definitions following Eq. (j2"7| we have 



irr(n) = (n, T) 
irr(T) = (n,T) 
irr(-u) = (u, v) 
irr(w) = (u, v) 



rev(n) = (it, v) 
rev(T) = (u, v) 
rev(w) = (n, T) 
rev('u) = (n, T), 



(34) 
(35) 
(36) 
(37) 



and the set of irreversible and reversible coefficients of the matrix are those with indexes given by, respectively: 

irr = (nn, nT, Tn, TT, uv, uu, vu, vv) (38) 
rev = (mi, nv, Tu, Tv, un, uT, vn, vT). (39) 

In the following we consider the two kinds of thermostat introduced above. 



1. Finite temperature thermostat, 75 > 
In the case "fb > 0, the dynamic matrix reads 

/ 



A = 



V 





n 






-(370W 





-Kk 2 






2 7 b) 






M\k 2 




7b) 









\ 



( 








—kn 












-Jfe2a 

na 







2 

n 


;„ Po 

mnT 








\ 











V 



(40) 



where v± and v\\ are the kinematic shear and longitudinal viscosity respectively, while k = 2k /(nd) is the thermal 
diffusion coefficient associated with the thermal conductivity k. Here 70 = (I — a 2 )/2d, gi is the pair correlation 
function at contact (in two dimensions we use the Verlet-Levesque approximation gi = (1 — 7p/lQ)/{l — p) 2 , po is 
the pressure, which for elastic hard disks reads po = nT(l + 2ng2nr 2 / A), and c is the thermal velocity. Here and 
in the following the transport coefficients v\\ and 7t are evaluated by the Enskog theory for elastic disks [33| at 
temperature T (which is the NESS temperature and may depend on inelasticity). Note that irreversible coefficients 
of A are always real, while reversible coefficients are always imaginary. 

Hydrodynamic noise is the sum of external noise due to the bath and internal one due to collisions. Only the latter 
is assumed to satisfy a "local equilibrium" assumption, i.e. the fluctuation-dissipation relation with the granular 
temperature [f9|, [2(|. Noises amplitudes have, then, the following amplitudes: 



D nn = 

D uu 





2k 2 KT 2 


dn 


nd 


lbT b ! v\ 


k 2 T 


n 


n 


lbT b ! v_ 


L k 2 T 


n 


n 



(41) 
(42) 

(43) 

(44) 



Following the general formula (|29[) . we get - for the entropy production rate at a given k: 



Ok 



Am i A TT A Tu i A uu A uT 

Dtt kn Dtt kn D uu 



SnT* + c.c. 



b.t. 



(45) 
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To obtain the above formula we have used the linearized continuity equation Sn = —i(kn)u which allows one to do 
the following replacements: u = iSn/kn and u* = —i5n jkn. Moreover, we have exploited the fact that all terms of 
the kind 4>4>* + c.c. are so-called "boundary terms" , i.e. they contribute to the entropy produced in the time-interval 
[0. t] by an amount |0(i)| 2 — </>(0)| 2 , which is usually sub-leading with respect to other terms which grow with time 
t [25l . [26l . [HI, [35]]. We have denoted all these terms by the symbol b.t. (those terms vanish in long-time averages, as 
dicussed below). Using the definition of the dynamical matrix A in Eq. (|45l) . we have 

a,, = h(k)R[Sn(k)f*(k)} + b.t., (46) 

where 

(T b - T) lb { Pq [k + 2(3/ d - lV||]fc 2 + g 2 n{ lb T b + is^T) + 6p Qlb T b /(dT)} 
( ) nT (Kk^T + 2 lb T b )(^T + lb T b ) ' ( V 

Here #t[x] denotes the real part of x and we have used the relation j uj ~ 2^ b (T b — T)/(dT) [36|. At equilibrium, 
namely for elastic interactions, T = T b and the entropy production vanishes identically, as expected. 

Notice that in the linearized hydrodynamics, modes at different k do not interact and therefore do not exchange 
energy or produce entropy. Entropy is instead produced separately at each k, because of unbalanced fluxes among 
different components of the fields, i.e. an entropy production rate exists for each mode k. The total rate of entropy 
produced by the hydrodynamic degrees of freedom is the sum of all exfc in the hydrodynamic range of ks. Interestingly 
enough, because of the linear approximation adopted, each scale k has its own entropy production which is expected 
to satisfy the Fluctuation Relation. 



2. Average entropy production 



Eq. (|46|) is the central result of this paper: it expresses the entropy production defined microscopically in terms of 
macroscopic quantities, following the structure of linear non-equilibrium thermodynamics. In particular, the quantity 
Sn(k)T*(k) plays the role of a current of energy, whereas the quantity h(k) is a thermodynamic force which is 
vanishing at equilibrium. Moreover, taking the stationary average over the stationary distribution on equation (|46l) . 
the boundary terms vanish, and, exploiting the dynamical equation for the temperature field, the entropy production 
can be expressed as a linear combination of static correlation functions (structure factors) 

K) = h(k)3t[(5n(k)t*(k))] = h(k){A* Tn C nn (k) + A* TT C nT (k) + A* Tu C nu (k)}, (48) 

where Cy-(fc) = {<fii(k)<j)*(k)) . The structure factors can be computed analytically. Then, for small k, we have 
h(k) ~ ao + a 2 k 2 , with 



(Gpo + dg 2 nT)^o /An ^ 



and $t[{6n(k)T*(k))] ~b 2 k 2 , so that 



(<7fc=o) = 0. (50) 

In the opposite limit, k — > 00, one finds the behavior (ufc) ~ fc~ 4 . In Fig. [TJwe show the average entropy production 
as a function of k for several different values of the parameters a and p. The first observation is that the average 
entropy production defined in Eq. (j48|) vanishes by definition at equilibrium. In this situation, the coefficient A^ n is 
zero and the correlation functions C„t(A) and C nu between fields of opposite parity under time reversal, vanish due 
to time reversibility. Moreover, from the results of our calculations we find that, even in inelastic cases, both at very 
large and very short length scales, the system looks like in equilibrium, where we have a zero entropy production. 
From the results of a previous work [20j | , we know that modes at k — > show equipartition at the bath temperature 
T b , whereas modes for k very large at the granular temperature T. It makes sense that at the length-scales where a 
good degree of equipartition is attained the entropy production turns to be zero. Another interesting similarity with 
the results on fluctuating hydrodynamics obtained in [2(| is that, also here, from the study of the entropy production, 
a characteristic non- equilibrium length-scale emerges. It seems that there is a preferential wave-vector at which the 
excitation of modes in the system produces the higher degree of dissipation, and that this wave- vector corresponds to 
a finite length. 
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FIG. 1. Average entropy production (tJk) as a function of kr, Eq. ((46} , for the driven granular model with finite temperature 
thermostat and friction yi > (defined by the Eqs. (|33|) and (|40I44|| ). The different curves show the behavior of {uk) for 
several values of the parameters a and p, as reported in the legend. Notice the peak at a certain wave-vector which signals the 
presence of a characteristic length-scale, as explained in the text. 



3. Infinite temperature thermostat 
For the case 7& = and T\, — > oo with finite T — I^T^, the dynamical matrix reads 



/ 



.4 







-7o 







and noise amplitudes are 







^f- -(3 7 o^ + ^ 2 ) 



\ ( 



1 



-kn 0\ 


mnT 



V o 



o 





1 /4Tr 4fc 2 «T : 



2 \ dn nd 

1 /T 2^||fc 2 T 

2 \n n 
T 2^,A; 2 T 



2 \n n 

Then, following exactly the same computation as before, we find 



0/ 



(51) 

(52) 
(53) 

(54) 

(55) 



(cr fe ) =h'(k)$t[(5n(k)f*(k))] 



with 



(56) 



h'(k) = - 



7o^o{[ 2 Po(« + ^||) + dg 2 nTv\\\k 2 + (6p Q + dg 2 nT)-i a uj a } 
2nT 2 (nk 2 + 2~f oj Q )(is\\k 2 + 7 w ) 



(57) 



In this case, the series expansion of the prefactor h'(k), for small k, gives h'(k) ^ a' Q + a' 2 k , with 
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FIG. 2. Average entropy production {at) as a function of kr, Eq. (|56[) , for the driven granular model with infinite temperature 
thermostat and without friction, Tb = oo, 7t = with 2jbTb = T finite, defined by the Eqs. (|33|) and (|5H55|) . The different 
curves show the behavior of (0*.) for several values of the parameters a and p, as reported in the legend. Notice that reaches 
a finite value for k — > 0. 



u AnT 2 
while (3l[6n(k)T*(k)]) ~ b' + b' 2 k 2 , with 



6ft, + 4»nT 



b> = nT 2 (6p + dg 2 nT)j Q ujQ 

6po + dg 2 npoT + 9dn 2 Tj i>\\ujo ' 



so that for the entropy production at k = we have the finite value 



/ k = (6p + dg 2 nT) 2 ~f uj 

[ak=0 > A(6p 2 + dg 2 np T + 9dn 2 T lo ^u;o)' 1 j 

This result is qualitatively different from that obtained in the case of a finite temperature thermostat. The finite value 
for k obtained in this case could be related to the scale-free correlations in the fluctuations of the hydrodynamic 
fields of this model [37[- In the opposite limit, k — > oo, we again find (ofc) ~ fc~ 4 . In Fig. [2] we show the average 
entropy production as a function of k for several different values of the parameters a and p. Notice here that there 
is always a finite entropy production, even at very short wave- vectors, where a good degree of equipartition between 
modes is never obtained. The result we find is consistent with the fact that it does not exist any finite length-scale at 
which we can coarse grain the local hydrodynamics fields in order to see an equilibrium-like behavior. 



B. Homogeneous cooling 



In the limit r — > 0, the above model becomes a pure cooling granular flui d ( no energy injection). If the initial state 
is spatially homogeneous, a regime exists where homogeneity is preserved [38|, the temperature decays following the 
Haff's law and any time-dependence is enslaved by the temperature. This regime, known as Homogeneous Cooling 
State (HCS) lasts until instabilities arise and homogeneity is broken. In the HCS, it is useful to introduce a new time- 
scale, defined by dr = uj(t)dt where uj(t) is the time-dependent collision frequency and, at the same time, rescaling 
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density, velocity and temperature fluctuations by n (which is constant), y/2T(t) and T(t) respectively. With such 
prescriptions one achieves a stationary representation of the HCS whose Fourier-transformed linearized hydrodynamics 
has been studied for example in [U 39]. In this case, the dynamical matrix reads 



.4 



( ° 

-7o52 








-( 70 + ^fc 2 ) 






7o 



and noise amplitudes are 



70 - v\_k 2 J 









-fcA 










kc 2 Xo 

















D„ 



Dtt — ~ 



K 



2 nduj 

1 J/||fc 2 

2 nw 
lu±_k 2 
2 nw 








(61) 

(62) 
(63) 

(64) 

(65) 



Here k = n/w, = vit/ui and if± = v±/uj are the rescaled transport coefficients, Ao is the mean free path and 
c 2 = l/(2T)(dp/dn) T . Note that the above matrices are obtained by replacing 7fc — > —70^ in the matrices for the 
driven case (and of course applying the other rescalings discussed above). Note also that the real part of some of the 
eigenvalues of the hydrodynamical matrix become positive for wave- vectors smaller than the cut-off value fc*: such an 
instability, which marks the end of the HCS, can be prevented by taking a system small enough. Considering values 
of fc < fc*, for the HCS we find the following expression of the entropy production 



(o- fe ) = h"(k)^St[6n(k)f*(k)] 



(66) 



with 



h"(k) = _^±1^M, (67 ) 
2k«t^|| 

In Fig. [3] we show the behavior of the average entropy production as a function of fc/fc*, for different values of a 
and p. Notice that, even approaching the value fc*, where instabilities arise, the entropy production remains finite. 
Also in this case, in the limit fc — > 00 we have the behavior (<7k) ~ fc~ 4 . 



V. CONCLUSIONS 



The coarse-grained dynamics of spatially extended systems can often be described by a set of coupled Langevin 
equations. For this class of stochastic systems, we have presented a general formula for the fluctuating entropy 
production, showing that it can be put in a form similar to that expected from macroscopic theories, such as linear 
non-equilibrium thermodynamics. As an application, we have studied in detail the fluctuating hydrodynamic equations 
describing the large scales fluctuations of different kinds of granular fluids. In particular, we have singled out that, 
in homogeneously driven granular fluids, the continuous flow of energy across the system produces some correlations 
(SnT) , usually absent at equilibrium. Namely, a density fluctuation in a certain point of space couples to the variations 
in time of the local temperature in the neighborhood. We identified the quantity SnT with the out-of-equilibrium 
fluctuating current of such a system. Notice that, due to the homogeneous driving mechanism, this current has no 
direction in space (only in time), differently from many examples of off-equilibrium systems. Moreover, we have shown 
that the wave- vector dependent entropy production, , is a quantity that provides a measure of how much dissipativc 
is the granular system when variables are coarse-grained on a length-scale A <~ 1/fc. 
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k/k 

FIG. 3. Average entropy production (at) as a function of k/k* , Eq. (|66[) . for the homogeneous cooling state, defined by the 
Eqs. (|33|l and (|6U65[) . The different curves show the behavior of (<Tfc) for several values of the parameters a and p, as reported 
in the legend. Notice that a k reaches a finite value for k — k* , where k* is the wave- vector of the instability in the system. 



Appendix A 

In this Appendix we show that the entropy production vanishes identically when equilibrium conditions are con- 
sidered. In the following sum over repeated indices is always meant. 

Our purpose here is to skccth the main idea, and this is best done - for the sake of compactness and without losing 
in generality - restricting to the case of all real fields. For such a case the entropy production is given by 

a'{t) = a(t) + lnp(>(0)) - lnp(0(i)) = f dsB k , lrr ( S )[<j> k { S ) ~ B Krev {s)] + lnp(0(O)) - lnj>(0(*)) 



o 



n 



= ^kkV")* [&(*)&(*) - h(0)M0)] ~ D kk J d S (A* rr )k3<PM rev hi(f>i + lnp(0(O)) - lnp(0(i)), 

(Al) 

where, following the definitions (|21I22|) , (A irr ) l:j = l/2(A y + e^eA and (A re %- = l/2(A l3 - CiAijCj). Using the 
equilibrium distribution 

p{4>)^ex^[-^4>k{C- x )ki4>i\, (A2) 

together with the equilibrium relation 

A irr C = -D => D- x A trr = -C~\ (A3) 

one immediately finds that, in Eq. (|A1I) . the first term on the rhs of the last line exactly cancels the quantity 
]np((f)(0)) — hip((f>(t)). The remaining term can be recast in the form 

D k k(A irr ) kj ct>AA rev )ki^ = -{C- l ) kj h{A r ™)ki<t>i = -(C-^kfaiA™)^ = -4>C~ x A r ™^ (A4) 
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where the first equality follows from the relation (|A3|) . while the property C = C T has been used in the second. Now, 
since at equilibrium 

A rev c + C ( A revjT = Q => fT 1 A reV + (C~ l A rev ) T = 0, (A5) 

one has that (j~ 1 A rev is antisymmetric, namely 

4C- 1 A rev cj) = 0. (A6) 

This implies a 1 = 0. 

For the full complex case, the demonstration follows in a similar way. 



Appendix B 

In this appendix we present some complementary considerations in terms of the Fokkcr-Planck equation associated 
with the Langevin equations (|TB)) . This allows us to discuss relations involving the Onsager coefficients. 

The Fokker-Planck equation obeyed by the probability density for a linear process defined by Eqs. (jT6]) and (|27|) 
reads 

T = -^|:^^*)]' ( B1 ) 

with 

j(4>,t) = -A4-Rf(<f>,t). (B2) 
Here j is the probability current vector, A the dynamical matrix, 

dlnpQM) 

Ji^ t > = ' ( B3 ) 

and i? is any matrix which satisfies R s = D with R s its symmetrized and D the matrix of noise correlations (which is 
always symmetric by definition). This arbitraricty on R will be discussed below. We anticipate that, at equilibrium, 
R coincides with the matrix of the Onsager coefficients L. 

The stationary solution of Eq. (|B1|) . p{<j>) = lim^oo p(cj>, t), satisfies 

hx\p((f>)] = const. - ]^<j)G4>\ (B4) 

where <jr is the adjoint of the vector (f> (transpose and conjugate) and G = C _1 with Cjj = (4>i4>*j) the self-adjoint 
matrix of covariances. Let us also define the quantity 

f {<!>)= lim f{^t) = -G4>. (B5) 
By direct substitution one finds that C must satisfy the following equation: 

In general this means that 

AC = R + Q (B7) 

with Q an antihermitian matrix, i.e. Q = —Q^. The relation with the D matrix is also ambiguous, i.e. AC = D + Q' , 
with Q' another (in principle different from Q) antihermitian matrix. In the stationary state one has, for the probability 
current: 

j°(0) = lim j(t) = —{RG)<j> - {QG)<j> + (RG)<j> = -{QG)<f>. (B8) 

t— VOO 



15 



Therefore, the probability current at any time may be rewritten as 



i?[/(<M) -/V)] +/(</>)■ 



(B9) 



The equilibrium condition is j (<f>) = 0, which requires Q = 0, leading to R = AC, which is the fluctuation- 
dissipation relation (FDR). Moreover, at equilibrium, from the time translation and time inversion invariance of the 
joint distribution p((j)(t),t; 0(0), 0) = p((p(0),t; 4>{t), 0), there follow the relations 



On the other hand, a non-equilibrium stationary state with j ^ can be observed for open systems where energy or 
matter is exchanged with the boundaries/thermostats with a preferred direction in time: this is the case of the driven 
granular fluids discussed in Section HVl 

Finally, note that in the Fokker-Planck equation, which is uniquely determined by the linear Langevin equation, 
only A and D = R s are fixed (and consequently also C or G). The matrix R is determined only in its hermitian 
part, while its antihermitian part is undetermined. At equilibrium this indeterminacy disappears, because of the FDR 
which fixes R = AC. Note also that the FDR does not imply R (or AC) to be symmetric. Its symmetry properties 
are determined by the parities under time-reversal (|B11I) . as discovered by Onsager [23|,[24|: Rij = Rj^ej. Moreover, 
let us remark that the Onsager matrix is uniquely defined only at equilibrium. Indeed one has at least two options 
to define the Onsager matrix in the non-equilibrium case: 1) one can follow the principle that the Onsager matrix is 
the "proportionality factor" between thermodynamic force f — f° and the current excess with respect to j°, so that 
Eq. (|B9[) implies L = R; 2) one can follow the principle that the Onsager matrix is the proportionality factor between 
the decay of fluctuations and the stationary entropic force X — —f°, i.e. d(<fi)/dt = —L(X), and this implies L = AC. 
At equilibrium the two options coincide. Of course those two principles are still arbitrary: Onsager matrix is, indeed, 
undefined in non-equilibrium dynamics. 
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